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ABSTRACT 


This paper presents a novel cylindrical-shaped Flexure-based Electromagnetic Linear Actuator (FELA) 
that exhibits predictable and re-configurable open-loop positioning resolution. By combining contactless 
Lorentz-force actuation and frictionless flexure-based supporting bearings, it produces high repeatable 
motion and sub-micron positioning resolution. In this paper, the design concept of this cylindrical-shaped 
FELA will be introduced. It focuses on the modeling of the flexure-based supporting bearings, the thermal 
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modeling of the electromagnetic module, and the unique characteristics of FELA, i.e., predictable and re- 
configurable open-loop positioning resolution. A prototype was developed to evaluate the performance 
and demonstrate these unique characteristics of this new class of nanopositioning actuator. 


1. Introduction 


An actuator with nano-positioning capabilities has been the 
key technology in nano-/micro-scale manufacturing processes 
such as nano-imprint lithography, fibre optics alignment, MEMS 
assembly, nano-scale machining, etc. Among various types of nano- 
positioning actuators, the solid-state piezoelectric (PZT) actuators 
have been the most popular choice due to their large actuating 
force and high stiffness. However, PZT actuators have limited dis- 
placements that make them unsuitable to drive high-precision 
manipulators targeted for large traveling range | 1]. Although some 
existing high-precision positioning actuators are able to eliminate 
such limitations, the displacement amplification techniques that 
are used within these actuators inherit other drawbacks. For exam- 
ple, PZT-driven actuators that use high-pitch screw actuating-shaft 
to achieve millimeters of displacement have poor repeatability due 
to backlash and Coulomb friction [2]. Others that use the mag- 
netostrictive clamping technique [3], the inchworm clamping [4], 
and the impact-force method [5] to drive an internal shaft for 
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displacement amplification purposes have low payload capacities. 
In addition, the slow response speed makes these actuators ineffi- 
cient for high speed applications. 

Electromagnetic driving scheme has the potential of deliver- 
ing millimeters of traveling range with nanometers of positioning 
resolution. Actuators of such frictionless drive, i.e., voice-coil lin- 
ear actuators and solenoid actuators, have also been employed in 
high-precision manipulation [6,7 ]. However, the voice-coil actuator 
produces relatively small output forces (or poor force sensitivity). 
Although the moving magnet actuator offers good dynamic behav- 
ior and good heat dissipation, the magnetic force attraction due 
to the external iron casing causes the output force to be incon- 
stant throughout the entire displacement stroke [8]. A solenoid 
actuator offers good force sensitivity but is unable to achieve a con- 
stant output force throughout its allowable traveling range due to 
the nonlinearity between force and displacement [9]. Furthermore, 
magnetization on the ferromagnetic stator introduces nonlinear 
magnetic hysteresis that leads to inaccuracy in nanometric posi- 
tioning tracking [10]. The drawbacks of existing nano-positioning 
actuators motivated the development of a Flexure-based Electro- 
magnetic Linear Actuator (FELA) [11]. FELA is formed through a 
marriage between an electromagnetic (EM) driving scheme and 
the flexure joints to achieve a few millimeters of displacement, 
large continuous thrust force, and a direct-force control capability. 
From past literatures [12,11], FELA was able to deliver a positioning 


accuracy of +20 nm, a continuous output force of 60 N/A, an actuat- 
ing speed of more than 100 mm/s throughout a large displacement 
stroke of 4mm. 

In our previous research efforts, the thermal management of 
FELA was not considered because its applications were mainly 
on automating embossing [13,14] and nano-imprinting processes 
[15,16] to replicate micro/nano-scale features. In these processes, 
large output force with direct precise force feedback is essen- 
tial. Hence, high force generation and high motion repeatability 
were the key considerations during the design phase. In addition, 
the closed-loop feedback control implemented on the end-effector 
ensures any form of thermal expansion will not affect the position- 
ing accuracy. Thus, the coils usually operate between 60 and 80°C 
during those imprinting operations. However, recent challenges 
arisen from applications such as the optical fiber alignment and the 
bio-cell manipulation have demanded FELA to produce higher posi- 
tioning and thermal stability over a few minutes of operation period 
with or without closed-loop feedback control. Thus, an improved 
version, which could meet such requirements, is required for those 
applications. In our recent research efforts, a model-based ther- 
mal compensation control was implemented to the existing FELA 
to compensate for the material expansion at the tip of the out- 
put shaft due to thermal induction from the energized coil [17]. 
This investigation shows that the effectiveness of thermal con- 
trol largely depends on the amount of thermocouple sensors used 
to estimate the thermally induced position error. To reduce the 
complexity of the entire control system, a two-stage optimization 
method was explored to re-design the Electromagnetic (EM) mod- 
ule of the existing FELA [18]. These efforts led to the conclusion 
that by factoring the thermal effect in the initial design stage will 
be much effective over the thermal control or thermal management 
approach in the later stages. 

This paper presents a novel cylindrical-shaped FELA that 
achieves lower heat generation as compared to the existing 
rectangular-shaped FELA. It consists of anew Lorentz-force Electro- 
magnetic (EM) module and flexure-based membrane supporting 
bearings. In this work, the cylindrical-shaped FELA has targeted 
specifications of +10nm positioning accuracy and +0.15°C ther- 
mal stability over a stroke length of 2mm at the end-effector. 
To achieve the desired thermal characteristic, an accurate ther- 
mal modeling of the EM module will be presented. To synthesize 
the desired stiffness characteristic that facilitate the EM module in 
achieving the targeted thermal stability, the stiffness modeling of 
the flexure-based membrane bearings was conducted using a novel 
semi-analytic modeling approach [19] and will also be presented 
in this paper. Most importantly, this paper also presents a unique 
characteristic of the FELA, i.e., predictable and re-configurable 
open-loop positioning resolution. Such characteristic cannot be 
found in existing nanopositioning actuators. All theoretical mod- 
eling and the unique characteristics of FELA will be evaluated and 
demonstrated. 


2. Design concept of the novel cylindrical-shaped FELA 


The cylindrical-shape FELA comprises of a new Lorentz-force EM 
module and a pair of flexure-based membrane support bearings 
as shown in Fig. 1. Termed as an Electromagnetic Driving Mod- 
ule (EDM), the Lorentz-force EM module is formed by a PM-based 
magnetic circuit with the moving air-core coil while while the 
flexure-based membrane bearings are used to support both ends 
of the moving air-core coil in order to retain the contactless nature 
of the EM driving scheme. Consequently, the frictionless character- 
istic of both driving and supporting elements ensure high motion 
repeatability. 
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Fig. 1. A detailed breakdown on the cylindrical-shaped FELA. 


2.1, A segmented dual-magnet configuration 


In the previous EDM design [11], 50% of the moving air-core coil 
does not operate within the magnetic field regions. Assuming that 
this portion of coil could operate within the magnetic field regions, 
the overall force generation could be enhanced by up to 50% with 
the same amount of input current. In other words, only 50% of the 
input current would be needed to generate the equivalent amount 
of force that the previous EDM design is producing. Consequently, 
heat generation could be reduced by half. In this paper, a seg- 
mented Dual-Magnet (DM) configuration was proposed to ensure 
that the entire air-core coil is operating within the magnetic field 
regions as shown in Fig. 2a. Instead of a complete stator casing, the 
new cylindrical-shaped EDM was constructed by a group of seg- 
ments whereby each segment was formed via a DM configuration 
as shown in Fig. 2b. 

Based on a segmented architecture, the magnetic field travels 
from outer PM to the inner PM within the designated closed-loop 
ferrous path to reduce magnetic leakages. (Note: the magnetiza- 
tion direction of the outer and inner PMs are similar.) Thus, most of 
the magnetic field could be extracted from both PMs to enhance 
the magnetic flux density within the effective air gap, in which 
the moving air-core coil operates. The segmented concept also pre- 
vents demagnetization between two inner PMs or two outer PMs 
when packed closely together because of the designated closed- 
loop ferrous path from each segment. Such a concept also reduces 
the assembly time significantly since the PMs can be glued in each 
segment concurrently before assembling all the segments together. 
Lastly, the cylindrical DM configuration forms an encasement for 
the magnetic field and prevents magnetic field leakage to the envi- 
ronment. Hence, FELA can be used in certain applications and 
environment that are sensitive to EM or magnetic field disturbance. 

Unlike conventional magnetic circuits [20,21] that deliver 
inconsistent and non-uniform magnetic flux density within the 
effective air gaps, the DM configuration delivers constant and 
evenly distributed magnetic flux density within a large effective 
air gap [22], i.e., 10 mm gap between the outer and inner PMs. By 
restricting the moving air-core coil to operate within the effective 
air gap as shown in Fig. 3, the new cylindrical-shaped EDM produces 
a constant current-force sensitivity (N/A) throughout the entire 
traveling range of the moving coil. Governed by the Lorentz-force 
principle and assuming that the magnetic flux density, Bext, is per- 
pendicular to the direction of input current, i, the output force, F, is 
expressed as 


F = iLBext (1) 


Segment 


4 Magnetization direction 


qe 


Ferrous path 


Outer PM 


(b) 


Fig. 2. (a) A cylindrical-shaped EDM constructed by a group of segmented DM configuration (US patent appl. #2013/990,219). 


where L represents the total coil length operating within the effec- 
tive air gap. Aided by the segmented DM configuration, Eq. (1) 
also suggest that the new cylindrical-shaped EDM provides a linear 
current-force relationship due to the constant magnetic flux den- 
sity within the effective air gap. The modeling and optimization 
of the DM configuration have been presented in recent literatures 
[18,22]. This paper mainly focuses on the thermal modeling of the 
EDM. 


2.2. Flexure-based membrane support bearings 


To retain the contactless nature of the Lorentz-force EM driving 
scheme, flexure joints are used to develop frictionless bearings to 
support the moving air-core coil as illustrated in Fig. 3. Unlike other 
frictionless bearings such as the air bearings, the magnetic bearings, 
and the hydrostatic/hydrodynamic bearings, these flexure-based 
bearings do not require any air/electrical/fluid source, expensive 
sensors, nor complex control systems to function. As a result, 
flexure-based bearings are simple, inexpensive and maintenance- 
free. In addition, a well-designed flexure-based bearing ensures low 
stiffness in the actuating direction while providing relatively higher 
stiffness in all non-actuating directions within the allowable travel 
range. Unlike the previous FELA design, which used clamping blocks 
and stainless steel shims to form such flexure-based bearings | 11], 
the new cylindrical-shaped FELA adopts a flexure-based membrane 
architecture that is formed by a single monolithic-cut stainless 
steel sheet. Consequently, the absence of the clamping blocks also 
eliminates these secondary moving masses. As the flexure-based 
bearings only operate within the elastic region of the material, the 
relationship between the force generated from the EDM and the 
output displacement, 6, of the air-core coil translator is assumed to 


be linear and expressed as 
F 
K=5 (2) 


where K represents the stiffness of the flexure-based bearings. 
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Fig. 3. Working principle of FELA: flexure-based bearings supporting a moving air- 
core coil that operates within the effective air gap of a DM configuration. 


2.3. Re-configurable open-loop positioning resolution 


Due to the linear current-force relationship of the EDM and 
the linear force-displacement relationship of the flexure-based 
bearings, the open-loop positioning resolution of FELA becomes 
predictable and re-configurable. (Note: an open-loop positioning 
resolution is defined as the achievable resolution based on the 
smallest controllable input current.) By combining the Ohm’s Law, 
i.e., V=/R, with Eqs. (1) and (2), the open-loop positioning resolution 
of FELA is given as 


$a V(Reoit + Ree) a (3) 


where R,oi represents the resistance of the air-core coil and Rext 
represents the resistance of the external resistor. Eq. (3) suggests 
two different approaches to enhance the open-loop positioning res- 
olution of FELA. First, by increasing the stiffness of flexure-based 
bearings to enhance the positioning resolution for the same amount 
of driving force. Second, by adding external resistors between the 
voltage supply and the system to reduce the amount of current 
via the voltage mode. Flexibility to reconfigure the open-loop posi- 
tioning resolution is an unique characteristic of FELA and will be 
demonstrated in later part of this paper. 


3. Thermal modeling of the cylindrical-shaped EDM 


Thermal modeling and analysis plays a crucial role in syn- 
thesizing the cylindrical-shaped EDM to achieve a desired power 
consumption and temperature characteristic for continuous force 
generation. In this work, the lumped capacitance modeling 
approach was used to obtain the mathematical representation of 
the thermal performance of the EDM. 


3.1. Power management 


From Eq. (1), the output force is a function of input current, mag- 
netic flux density and the total coil length operating within the air 
gap. Based on the architecture of the cylindrical-shaped EDM, the 
coils are wound on a circular bobbin to form the moving air-core 
coil as shown in Fig. 4. Assuming that the coiling is ideal, the total 
number of rows will represent the total number of layers while 
each column will represent the number of turns per layer. Based 
on this configuration, the total coil length is represented as 


L = 27Nppi[rpN; + OD(N; + 1)(Ni/2)] (4) 
where N; represents the number of layers of coil within a bob- 


bin, Np represents the number of turns per layer, OD represents 
the overall diameter of the conducting coil wire (with insulation 


Z column Z, 
EE i 


Ideal coiling configuration 


Fig. 4. A circular bobbin that holds the moving air-core coil. 


layer), and rp represents the radius of the bobbin. Subsequently, 
the internal coil resistance of the total coil length is given as 


xL 
Rooi n(ID/2)° (5) 
where x represents the resistivity of the conducting material, i.e., 
1.7 x 10-8 Qm for copper wire, and ID represents the internal diam- 
eter or the core diameter of the conducting wire (excluding the 
insulation layer). With R,;, and input current, the power, Q, con- 
sumed by the EDM is expressed as 


Q = i? Reoil (6) 
3.2. Lumped capacitance modeling approach 


In nanopositioning applications, the air-core coil is usually slow 
moving and needs to hold at a position during operations. When 
the air-core coil is stationary, i.e., no work done, it can be assumed 
that all the power generated from the coil turns into heat. Thus, 
the elevation and saturation of coil temperature will reflect the 
thermal characteristic of the EDM, which is crucial for understand- 
ing the thermal stability of FELA. Based on the architecture of the 
cylindrical-shaped FELA, the moving air-core is located in the center 
while the stator and other structures constitute the thermal paths 
to dissipate heat to the external environment. The thermal perfor- 
mance can be represented as a classical analogous lumped-element 
circuit model that comprises of thermal resistors and capacitors 
[23]. Based on this approach, the thermal resistance of an element, 
R, is given as 

L 
err (7) 
where L represents the length of the thermal path, A represents 
the cross-sectional area of the thermal path, and h represents the 
thermal conductivity of the thermal path. In addition, the thermal 
capacitance of an element, C, is given as 


C= MH, (8) 


where M and H; represent the mass and the specific heat capacity 
of the element respectively. By differentiating the transient behav- 
ior of a resistance—capacitance (RC) circuit, the rate of change in 
temperature is expressed as 


AT=QR E exp (=)| (9) 


where t represents the time and t represents the thermal time 
constant, i.e., T = MH;R. 


3.3. Thermal model of the EDM 


Based on the lumped capacitance modeling approach, the ther- 
mal path of the cylindrical-shaped EDM is represented by two 
cascaded RC circuits as shown in Fig. 5. The source of this circuit 
is the temperature of the air-core coil, which is represented by Toi, 
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Fig. 5. RC circuit representation of the cylindrical-shaped EDM. 


and the ground represents the ambient temperature outside of the 
EDM, i.e., Tampient- The first RC circuit represents the thermal path of 
the air gap between the surface of the air-core coil and the surface 
of the outer PM where Rg represents the thermal resistance of the 
air gap while Cg represents the thermal capacitance of the coil. In 
addition, it also contains the parallel thermal path of heat transfer 
from the air-core coil to the bobbin that holds it where Rg repre- 
sents the thermal resistance of the air gap while Cg represents the 
thermal capacitance of the coil. For the first RC circuit, the thermal 
resistance and capacitance are expressed as 


Ri =(1/R¢+1/Rg) (10) 
Ci = Ce + Cg (11) 


Based on the first RC circuit, the change of temperature, AT}, 
due to the power generated by the input current is given as 


AT; = QR (12) 


The second RC circuit represents the parallel thermal paths from 
the outer PM to the environment and from the inner air gap to 
the environment. From the outer PM to environment, Rpyo and 
Rsop represent the thermal resistance of the outer PM and the sta- 
tor outer casing respectively while Cpyo and Csop represent the 
thermal capacitance of the outer PM and the stator outer casing 
respectively. From the inner air gap to environment, Rc, Rpyy, Rsp, 
and Rss represent the thermal resistance of the inner air gap, inner 
PM, stator inner shaft, and stator side structure respectively Simi- 
larly, C represents the thermal capacitance of these elements. For 
the second RC circuit, the thermal resistance and capacitance are 
expressed as 


1 
Ry = [1/(Rpmo + Rsop) + 1/(RG + Remi + Rsip + Rss)] + Rconv 


(13) 


Cy = (1/Cpmo + 1/Csop)* + (1/Cpmi + 1/Csp + 1/Css) | (14) 


Based on the second RC circuit, the change of temperature, AT>, 
due to the power generated by the input current is given as 


AT) = QRy (15) 


A generic temperature-time characteristic of the air-core coil 
due to a constant current input over a period of time is plotted 
in Fig. 6. AT,, which is contributed by the first RC circuit, occurs 
rapidly because the thermal path is short due to the small air gap 
between the coil and the PM. In addition, the rise is linear as the 
thermal capacitance of the bobbin is low due to its small volume 
mass that is unable to store much heat. Hence, the heat generated 
from the coil is quickly dissipated to the surrounding PMs. In most 
cases, AT, would be at 95% of its final value by 3t. On the other 
hand, the temperature rise from AT, to AT> is usually longer and 
nonlinear as compared to AT}. This is because there are more ele- 
ments in the second RC circuit that lead to higher thermal resistance 
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Fig. 6. A generic temperature-time graph of an air-core coil when introduced with 
a constant current over a period of time. 


Table 1 
Selected values for the thermal properties of different elements. 
Elements Thermal properties Values 
Copper (Coil Density (kg/m? 8940 
wire) Specific heat capacity (J/kg K) 390 
Thermal conductivity (W/m K) 400 
Iron (Stator) Density (kg/m? 7850 
Specific heat capacity (J/kg K) 449 
Thermal conductivity (W/m K) 80 
Rare-earth Density (kg/m? 7400 
(PMs) Specific heat capacity (J/Kg K) 506 
Thermal conductivity (W/mK) 17 
Aluminum Density (kg/m? 2700 
(Non-ferrous Specific heat capacity (J/Kg K) 900 
parts) Thermal conductivity (W/m K) 250 


and capacitance. Higher thermal resistance will cause the temper- 
ature rise to be slow and higher thermal capacitance will cause 
nonlinearity. In most cases, it will take more than 3t before AT 
reaches 95% of its final value. Without considering the transient 
temperature over a period of time, the final saturation temperature 
of the air-core coil is expressed as 


Thinal = Tambient + AT; a AT? (16) 


Lastly, Table 1 lists the selected values for the thermal proper- 
ties of different elements that are used to calculate the saturation 
temperature of the air-core coil based on the thermal resistance 
and capacitance of the proposed RC circuit. 


4. Stiffness modeling of the flexure-based membrane 
support bearings 


After designing the cylindrical-shaped EDM based on the pre- 
sented thermal modeling, the predicted output force at the desired 
operating temperature will be used to design the flexure-based 
membranes. Fig. 7 illustrates a flexure-based membrane support 
articulated by four parallel beam-based flexure joints. The length 
of each beam-based flexure joint is a combination of four folds of 
flexible beams, which are coupled together to deliver a pure trans- 
lational motion. With the fixed segment attached to the stator of 
the EDM, a single-DOF motion is produced at the moving segment 
through the deflection of these beam-based flexure joints. 


4.1. Semi-analytic model 


The stiffness of the flexure-based membrane can be obtained 
using the semi-analytic model [19]. This theoretical model provides 


a simple and generic solution for approximating the large non- 
linear deflection of any beam-based flexure configuration, i.e., a 
beam-based flexure coupled with a rigid-link of any length, shown 
in Fig. 7d. The accuracy and robustness of the proposed model 
is ensured by considering two crucial factors while adopting the 
principle of using a torsional spring to represent the stiffness of a 
beam-based flexure joint. The first factor is the parasitic shifting of 
the ‘pivot’ joint of a beam-based flexure joint during large deflec- 
tion, while the second factor is the changing angular stiffness of the 
torsional spring. With all derivations being presented in previous 
work [19], the semi-analytic model states that the deflection of a 
flexure configuration is given as 


l 
ô= | L+ —)si 17 
( + z5) sing (17) 
where « represents the deflection angle and @ is a Sinc function, 


i.e., W =sina/a. The resultant parasitic motion, ôp, normal to the 
deflection axis is given as 


bp = ¢ H 5) ¢ H oa) cosa (18) 
The loading force, F, is expressed as 
Ela (19) 


i (L+ (pl/2@)) sin((2r/2) — a) 


where E and I represent the Young’s Modulus and second moment 
of area of the flexure joint respectively, while p is 


1V/1.8+L 
= ee) 
The maximum bending stress, Omax, is given as 
F[l/2+(L+1/2@)cosa](h/2 
ig UE eae as KI2) ei) 


where h represents the thickness of the flexure joint. 

(Note: Termed as semi-analytic model since deflection or load- 
ing force can be solved directly with a given deflection angle 
while solving deflection angle based on a known deflection or 
loading force would require iterative root-finding procedures, e.g., 
Newton-Raphson method, etc.) 


4.2. Modeling of the flexure-based membrane 


In this work, the flexure-based membranes are made up of 
beam-based flexure joints with no rigid-link coupled to them. 
Hence, L=0 when using the semi-analytic model. The deflection 
stiffness along the x-axis of a beam-based flexure joint with no 
rigid-linkis denote as Ks. Hence, the loading of the beam-based flex- 
ure joint becomes a point force. However, the point force loading 
does not permit a pure translational motion. As each flexure mem- 
brane is a supporting bearing for the air-core coil translating in a 
single linear motion, the beam-based flexure joint must be modeled 
as a cantilever beam with a pure translation motion. To deliver a 
pure translation motion, each beam-based flexure joint deflects in 
”S”-like shape and can be represented as two identical beam-joints 
with individual length being 1/2 and the desired translation motion, 
A, becomes twice of the deflection of each beam-joint, i.e., A/2, as 
shown in Fig. 7c. Hence, the translation stiffness is given as 


Kj = 25 (22) 


where driving force, F,, is obtained from Eq. (19) using the deflec- 
tion angle, a, derived from Eq. (17) based on ô= A/2, | 1/2, and 
L=0. With each limb being formed by two identical beam-joints 
connected in series, the translation stiffness along the x-axis of 
each limb, K;, is half of a beam-joint, i.e., K; = K;/2. Consequently, 
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Fig. 7. (a) A flexure-based membrane support bearing articulated by four parallel beam-based flexure joints. (b) Each beam-based flexure joint combines four folds of flexible 
beams to deliver (c) a pure translational motion. (d) Illustration of a generic beam-based flexure configuration. 


the translation stiffness of a pair of flexure-based membrane is 
governed by 


Ks = 8K, (23) 
As mentioned, the length of each beam-based flexure joint is 


a combination of four folds of flexible beams. The length of each 
flexible beam is represented as 


l = (90° —d) (=) r (24) 


where j represents the beam number {i= 1, 2, 3, 4}, dis the angle (in 
degree) and r is the radius of each arc that represents the length of 
the beam as shown in Fig. 7b. From the fixed segment to the mov- 
ing segment, d=2(d3) and r=r1 to determine the length of flexible 
beam #1, d=d2+d3 and r=r2 to determine the length of flexible 
beam #2, d=2(d2) and r=r3 to determine the length of flexible 
beam #3 and d=d1+d2 and r=r4 to determine the length of flexi- 
ble beam #4. Subsequently, the length of each beam-based flexure 
joint is given as 

4 
l= Soi (25) 

j=l 


J 


5. Prototype 


The presented analytical models were used to synthesize the 
cylindrical-shaped FELA based on the targeted specifications, i.e., 
+20nm positioning stability and +0.15°C thermal stability over a 
stroke length of 2mm at the end-effector. In addition, the overall 
dimension must not exceed @50 x 100 mm? due to size constraint 
of the entire positioning system. Based on the magnetic field 
model derived in [22,18], the thickness of the PMs within the seg- 
mented DM configuration were selected as 6mm with a length of 
50 mm. Using the presented thermal model, the cylindrical-shaped 
EDM was predicted to produce a constant current-force sensitiv- 
ity of 10.15 N/A. To achieve a thermal stability of +0.15°C, the 
stiffness of a pair of supporting flexure-based membranes plays 
a crucial role as the EDM must operate near ambient tempera- 
ture, i.e., 20°C. Using the presented model to synthesize the pair 
of flexure-based membrane support bearings, the bearings were 
eventually designed with a predicted stiffness of 1375N/m. As a 
result, only 0.136 Ais required to generate 1.38 N to achieve +1 mm. 
Based on the thermal model, the temperature of the air-core coil 
within the cylindrical-shaped EDM will rise by 0.7°C after 1.5h. 
As this predicted performance can satisfy the targeted specifica- 
tions, acylindrical-shaped FELA prototype based on the synthesized 
EDM and flexure-based membranes was developed as shown in 
Fig. 8. Table 2 lists the parameters and actual specifications of the 


Fig. 8. Prototype of the new cylindrical-shaped FELA. 


Table 2 
Parameters and specifications of the cylindrical-shaped EDM and the flexure-based 
membrane support bearing. 


Specifications 


Cylindrical-shaped EDM 
Wire diameter, OD 
Coil resistance, Reoil 


0.48 mm (insulation layer ~ 0.03 mm) 
5.93Q 


Flexure-based membrane 


Material SUS301 Full-hard 
Young’s modulus 193 GPa 
Thickness 0.2mm 

Width 3mm 


d1=35°, d2=25°, d3=20° 
r1 = 18.20 mm, r2= 14.90 mm, 
r3 = 11.60 mm, r4=8.30mm 


Angle of arc, d 
Radius of arc, r 


cylindrical-shaped EDM and the flexure-based membrane support 
bearing. 


6. Experiments and results 
6.1. Thermal characteristic of the cylindrical-shaped EDM 


An experiment was conducted to investigate the thermal char- 
acteristic of the cylindrical-shaped EDM and to evaluate the 
accuracy of the proposed model. A thermocouple sensor (RT-100) 
was attached onto the surface of the air-core coil and connected 
to a computer to log the temperature at an interval of 1s. Tem- 
perature measurement was taken only with the cylindrical-shaped 
EDM because the thermal characteristic was modeled without the 
flexure-based membranes and supporting components. To prevent 
heat loss through any material that is not part of the prototype, it 
was placed on a carbon fiber composite board, which was treated 
as an isolator. A set of constant current inputs was used to ener- 
gize the air-core coil and temperature of the coil was logged until 
it saturated, i.e., about 1.5h after being energized. The predicted 
and measured coil temperature based on the constant current val- 
ues are listed in Table 3. A comparison between the predicted and 
measured temperature values shows that the presented analyt- 
ical model is accurate in predicting the coil temperature with a 


Table 3 
Predicted and measured coil temperature after 1.5 h based on constant input current. 


Current (A) Power (W) Coil temperature after 1.5 h (°C) 
Predicted Measured 

0.1 0.059 20.4 20.3 

0.4 0.949 26.2 28.0 

0.5 1.483 29.6 32.0 
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Fig. 9. Current-force relationship of the cylindrical-shaped EDM due to different 
positions of the air-core coil translator. 


maximum deviation of 7.5% obtained the investigation. The slight 
deviation is contributed by variation in the internal coil resistance 
where actual value is measured at 5.93 Q while the predicted value 
is 5.31 Q. In addition, the selected values for each thermal property 
listed in Table 1 may differ from the actual values. Nevertheless, the 
experimental results have shown that the developed cylindrical- 
shaped EDM has achieved the desired thermal characteristic. 


6.2. Current-force sensitivity of the cylindrical-shaped EDM 


To evaluate the current-force sensitivity, the FELA prototype 
was mounted on a linear guide that locked the cylindrical-shaped 
EDM at any position along the displacement range. A force sensor 
(ATI, Nano40) was attached to the output shaft with the other end 
of the sensor mounted on a fixed wall. The flexure-based bearings 
ensure that no friction force will be contributing the force mea- 
surement. When the air-core coil was energized by an amount of 
current, the output shaft moved against the force sensor, which 
picked up the amount of force generated. The linear guide with a 
locking mechanism allows the cylindrical-shaped EDM to change 
its position with respect to the fixed air-core coil. As a result, the 
current-force sensitivity at different positions can be evaluated. 

The selected positions along the displacement range to evalu- 
ate the current-force sensitivity were neutral position, +0.5 and 
+1 mm. At each position, a power supply source was used to ener- 
gized the air-core coil from 0 to 1 A with an interval of 0.1 A. Fig. 9 
plots the input current and output force when the air-core coil was 
at all selected positions w.r.t. the fixed EDM. It shows that the rela- 
tionship between input current and output force obtained from 
all five positions are very consistent and the deviations between 
the experimental and predicted results are less than 8.35%. Table 4 
summarizes the current-force sensitivity at all selected positions. 
Results show that the current-force sensitivity at various positions 
along the displacement range are consistent with a deviation of 
0.58% recorded between the maximum and minimum values. Sub- 
sequently, the average current-force sensitivity throughout the 
entire stroke measured from the developed prototype is 10.3 N/A. 
The achieved current-force sensitivity is close to the predicted 
value of 10.2 N/A. Most importantly, this investigation shows that 


Table 4 
Current-force sensitivity at various positions along the traveling range. 


Pos (mm) —1.0 —0.5 Neutral +0.5 +1.0 
Current-—force (N/A) 10.36 10.34 10.30 10.32 10.36 
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Fig. 10. Experimental results and theoretical prediction on the stiffness of a pair of 
flexure-based membranes. 


the developed cylindrical-shaped FELA has achieved a near con- 
stant current-force sensitivity throughout 2 mm stroke. 


6.3. Stiffness of the flexure-based membrane support bearing 


Stiffness of the flexure-based membrane support bearing gov- 
erns the relationship between the force generated from the EDM 
and the displacement of the output shaft. To evaluate the stiff- 
ness of these bearings, a high resolution current source (KEITHLEY 
224/2243) was used to energize the EDM and the generated dis- 
placement of the output shaft was measured by a high resolution 
optical linear encoder (MicroE Systems, M3500). In this work, high 
resolution current source was used because it provides very sta- 
ble current output. Hence, very stable step displacements can be 
obtained to reduce measurement errors. Subsequently, the force 
generated due to the input current was determined based on 
the current-force sensitivity of 10.3 N/A. With the measured dis- 
placement and predicted force, the actual stiffness of the pair of 
flexure-based membrane support bearings was obtained and plot- 
ted in Fig. 10. The predicted stiffness of 1375 N/m obtained from 
the semi-analytic model was also plotted against the measured 
values. Results show that the pair of flexure-based membrane sup- 
port bearings has an average stiffness of ~1398 N/m. Although 
the input current was restricted to 80mA due to the limitation 
of the current source, the trend suggests that the stiffness of the 
flexure-based membrane support bearings is linear throughout the 
targeted stroke range. Most importantly, this investigation shows 
that the semi-analytic model is accurate in predicting the stiffness 
with a maximum deviation of 7.15%. 


6.4. Reconfiguring the open-loop positioning resolution 


This section presents how the open-loop positioning resolu- 
tion of FELA can be reconfigured by changing the thickness of the 
flexure-based bearings or by adding an external resistor. 


6.4.1. Changing the thickness of the flexure-based bearings 

Governed by Eq. (3), a change in the stiffness of the flexure-based 
bearing will alter the open-loop positioning resolution. In Eq. (19), 
the second moment of area is expressed as 


_ bt? 
~ 42 


where b and t represent the width and the thickness of the flexure 
respectively. Based on Eq. (26), changing the thickness is the most 
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Fig. 11. Open-loop step displacements of FELA using flexure-based bearing with 
0.2 mm thickness. 


effectively approach in altering the stiffness of the flexure-based 
bearing. In this work, an investigation, which involved the FELA 
prototype, was conducted to demonstrate such flexibility. Based 
on the initial thickness of 0.2 mm, a set of input current, i.e., from 
i = 22 to 110A at an interval of 22 pA, via the high resolution 
current source, was used to energize the moving air-core coil so 
as to step the output shaft forward and return to neutral position. 
The step displacements of the output shaft were recorded by the 
high resolution optical encoder (5 nm/count) and plotted in Fig. 11. 
At 22 pA, a displacement of 172 nm was achieved and 345 nm was 
achieved at 44 uA. Hence, the open-loop positioning sensitivity was 
7.82 nm/pA. 

Next, the thickness of each flexure-based bearing was changed 
from 0.2 to 0.5mm at an interval of 0.1mm. For each thickness 
change, similar set of input current was used to drive the output 
shaft forward and return to neutral position. Based on a pair of 
flexure-based bearings with 0.3mm thickness, the semi-analytic 
model predicted a stiffness of 4641 N/m. For 0.4mm thickness, the 
predicted stiffness for the bearings is 11003 N/m. For 0.5 mm thick- 
ness, the predicted stiffness is 21485 N/m. Using the measured 
current-force sensitivity, i.e., 10.3 N/A, the step displacements of 
FELA were predicted and listed in Table 5 based on different bear- 
ing thickness and input current values. The experimental results 


Table 5 
Reconfiguring the open-loop step displacement of FELA using flexure-based bear- 
ings with different thickness. 


Current (WA) Predicted Measured (avg) Error (%) 
0.3 mm thickness: Step displacement (nm) 
22 45.43 42.69 6.42 
44 89.46 86.00 4.02 
66 133.83 130.69 2.41 
88 177.16 174.37 1.60 
110 221.18 212.43 4.12 
0.4mm thickness: Step displacement (nm) 
22 19.15 20.56 6.84 
44 37.59 37.66 0.17 
66 56.17 54.36 3.32 
88 74.74 74,98 0.31 
110 93.17 94.89 1.82 
0.5mm thickness: Step displacement (nm) 
22 9.88 10.00 1.17 
44 19.48 20.00 2.62 
66 28.92 30.00 3.61 
88 38.42 40.00 3.95 
110 47.86 50.00 4.28 
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Fig. 12. Open-loop step displacements of FELA supported by flexure-based bearings 
with 0.3 mm thickness. 
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Fig. 13. Open-loop step displacements of FELA supported by flexure-based bearings 
with 0.4mm thickness. 


show that the measured step displacements agree with the pre- 
dicted values. With a maximum deviation of 6.84% found between 
the measured and predicted values, this experimental investiga- 
tion shows that the semi-analytic model is accurate in predicting 
the stiffness for such kind of flexure joint configuration, i.e., L=0. 
Most importantly, results also show that the open-loop positioning 
resolution of FELA is predictable. 

Figs. 12-14 plot the actual step displacements of FELA sup- 
ported by flexure-based bearings with thickness of 0.3, 0.4, and 
0.5mm respectively. Referring to Fig. 12, the open-loop position- 
ing sensitivity of FELA was enhanced from 7.82 to 1.94nm/pA 
when the thickness of the flexure-based bearings changed from 
0.2 to 0.3 mm. With thickness of 0.4 mm, the open-loop positioning 
sensitivity was further improved to 0.93 nm//WA. Lastly, FELA pro- 
duced an open-loop positioning sensitivity of 0.45 nm/wA with the 
flexure-based bearings having a thickness of 0.5 mm. Experimen- 
tal results show that the open-loop positioning resolution of FELA 
is reconfigurable by changing the thickness of the flexure-based 
bearings. 

On the other hand, the achievable displacements were also 
affected by such changes. By comparing the final displacements 
plotted in both Fig. 12 and 14, a FELA, which was supported 
by flexure-based bearings with thickness of 0.3mm, achieved 
212.43nm at 110A. With the same amount of input current, 
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Fig. 14. Open-loop step displacements of FELA supported by flexure-based bearings 
with 0.5 mm thickness. 


the FELA, which was supported by flexure-based bearings with 
thickness of 0.5mm, achieved just 50 nm. This comparison sug- 
gests that although the open-loop positioning sensitivity of FELA 
can be enhanced by increasing the thickness of the flexure-based 
bearings, this approach requires more force to overcome the 
increasing stiffness of those bearings. Consequently, more heat 
will be generated due to the increase of input current. In addi- 
tion, changing the stiffness has a direct effect on the dynamics of 
the FELA. Assuming the moving mass remains unchanged, reduc- 
ing the stiffness so as to minimize heat generation will reduce 
the dynamic performance, i.e., natural resonant frequency. On 
the other hand, both dynamic behavior and open-loop position- 
ing resolution can be enhanced by increasing the stiffness of the 
flexure-based bearings. Most importantly, the results obtained 
from the experimental investigation have demonstrated that the 
open-loop positioning resolution of FELA is predictable and recon- 
figurable by changing the thickness of the flexure-based support 
bearings. 


6.4.2. Adding an external resistor 

An alternate approach to enhance the open-loop positioning 
resolution of FELA without affecting the heat generation is by mak- 
ing use the Ohm’s Law. From Eq. (3), adding an external resistor, 
Rext in series with FELA will reduce the amount of current energiz- 
ing the EDM via the voltage mode. Consequently, less amount of 
current will lead to less force generation and thus produce smaller 
step displacement. An experimental investigation was conducted 
to demonstrate how external resistance loads can enhance the 
open-loop positioning resolution. 

At initial conditions, i.e., flexure-based bearings with 0.2mm 
thickness, internal coil resistance, Roi) = 5.93 Q, and Rext = 0 Q, a set 
of input voltage, i.e., from VDC = 0.064 to 0.256 mV at an interval of 
0.064 mV was used to drive the output shaft forward and return 
to neutral position. The step displacements of the output shaft 
were recorded by the same optical encoder as plotted in Fig. 15. At 
0.064 mV, a step displacement of 172.87 nm was achieved by the 
FELA. At 0.256 mV, a step displacement of 345 nm was produced. 
Hence, the open-loop positioning sensitivity was 1347.66 nm/mvV. 
Next, Rext was increased to 5, 15, and 25 Q respectively. For every 
increment, the same set of input voltage was used to drive the 
output shaft while the optical encoder recorded the step displace- 
ments. Table 6 lists the measured step displacements of FELA due 
to different external resistance loads. 

Using Eq. (3), the current-force sensitivity of 10.3N/A, 
Reo =5.93 Q and stiffness of 1375N/m, the predicted step 
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Fig. 15. Open-loop step displacements of FELA at initial conditions. 
Table 6 


Reconfiguring the open-loop step displacement of FELA by adding different external 
variable resistance. 


Voltage (mV) Predicted Measured (avg) Error (%) 
Rex = 5 Q: Step displacement (nm) 

0.064 43.43 43.02 0.95 
0.128 88.37 89.23 0.98 
0.192 125.81 128.84 2.41 
0.256 168.50 171.99 2.08 
Rex = 15 Q: Step displacement (nm) 

0.064 23.22 25.11 8.17 
0.128 45.68 46.78 2.42 
0.191 68.90 70.94 2.96 
0.256 92.11 95.52 3.70 
Rext =25 Q: Step displacement (nm) 

0.064 15.73 15.52 1.29 
0.128 30.70 28.46 7.32 
0.192 46.80 44.77 4.36 
0.256 62.53 61.43 1:77 


displacements of FELA at different external resistance loads are 
also listed in Table 6. For every change of external resistance 
load, the experimental results show that the step displacements 
of FELA were well predicted by Eq. (3). The measured step dis- 
placements agree with the predicted values and the maximum 
deviation between them is 8.17%. Using the same set of input 
voltage, the experiment also show that the step displacements of 
FELA reduced when the external resistance load increased. Results 
have proven that the assumption of increasing resistance load to 
lower the input current, which in turn reduce the step displace- 
ment, is valid. When Rext = 5 Q, the measured open-loop positioning 
sensitivity became 362.81 nm/mV. When Rext=25 Q, the mea- 
sured open-loop positioning sensitivity became 245.78 nm/mV. 
Results show that the open-loop positioning sensitivity of FELA was 
increased from the initial 1347.66 to 245.78 nm/mV just by adding 
a 25 Q external resistor. Consequently, this experimental investi- 
gation has demonstrated that the open-loop positioning resolution 
is predictable and reconfigurable through an external resistance 
load. 

As compared to changing the thickness of the flexure-based 
bearings, there are two main advantages of using an external resis- 
tance load to configure the open-loop positioning resolution. First, 
this approach will not increase the force required to achieve the 
targeted displacement since the stiffness of the bearings remains 
unchanged. As a result, targeted displacement can still be achieved 
without increasing the input current. Thus, the thermal charac- 
teristics of the FELA remains unchanged. Second, the flexibility 
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Fig. 16. Open-loop step displacements of FELA when Rex =5 Q. 
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Fig. 17. Open-loop step displacements of FELA when Rex = 15 Q. 


is higher as compared to changing the thickness of the bearings. 
In real applications, adding an external resistance load in series 
with FELA is more flexible than changing the thickness of the sup- 
port bearings. In addition, thickness of these stainless steel shims 
are limited to certain range or specific thickness. The limitation 
of adding external resistance load is that this approach cannot 
reduce the noise level caused by the environmental disturbance. 
Figs. 16-18 plot the actual open-loop step displacements of FELA 
when Rext = 5, 15 and 25 Q respectively. At different external resis- 
tance loads, the average positioning noise is always +15 nm. From 
Fig. 12-14, the average positioning noise drops from +15 to +10 nm 
due to change in the stiffness of the flexure-based bearings. As a 
result, Fig. 14 shows more stable open-loop step displacements 
plot as compared to Fig. 18. Another limitation of adding exter- 
nal resistance load is that the overall power consumption of the 
system will increase. Although this increase in power consump- 
tion has no effect in the heat generation since the external load is 
outside of FELA, the wattage of the power source, which drives the 
entire system, needs to be increased accordingly to the added resis- 
tance. Although each approach has its limitations, both approaches 
are able to configure the open-loop positioning resolution of FELA. 
Choosing between two approaches will depend very much on the 
targeted applications and a good understanding of the performance 
trade-offs. 
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Fig. 18. Open-loop step displacements of FELA when Rex: = 25 Q. 
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Fig. 19. Positioning stability at 1 mm with PID servo-control. 


6.5. Evaluating the positioning and thermal stability 


Using the same optical encoder as the position feedback sen- 
sor of the output shaft, a Proportional-Integral-Derivative (PID) 
control was implemented to control the positioning of the FELA. 
Through the PID control scheme, FELA achieved a closed-loop posi- 
tioning resolution of +10nm at +1mm. A laser interferometer 
system (SIOS, resolution: 1nm/count) was used to evaluate the 
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Fig. 20. Operating temperature of coil at 1 mm. 


actual positioning accuracy at the tip of the output shaft. At 1 mm, 
Fig. 19 shows that the tip of the output shaft is capable of deliv- 
ering an average positioning stability of +10 nm. Fig. 20 plots the 
temperature of the coil when FELA was operating at 1 mm. It shows 
that the coil temperature remained at 20.25 °C with a thermal sta- 
bility of +0.15°C. Based on the measured flexure-based bearing 
stiffness of 1398 N/m and current-force sensitivity of 10.3 N/A, the 
estimated input current is approximately 0.135A (0.8 W). Based 
on this amount of input current, the thermal model predicted that 
the coil temperature will increase to 20.4°C by approximately 92 s 
(1.5 min). The predicted temperature and time almost agree with 
the coil temperature profile plotted in Fig. 20 with a deviation of 
0.15°C. 


7. Conclusion 


This paper has presented a novel nanopositioning actuator 
that achieves a large stroke range of +1mm with a position- 
ing stability of +10nm and a thermal stability of +0.15°C. To 
achieve such desired thermal characteristics, a thermal model of 
the Lorentz-force EDM based on lumped capacitance modeling 
approach was presented. Experimental investigations have shown 
that the presented thermal model is accurate in predicting the 
coil temperature. This paper also presented a novel semi-analytic 
modeling approach that was used to synthesize the proposed 
flexure-based membrane support bearings. Experimental eval- 
uations have shown that the proposed semi-analytic model is 
accurate in predicting the stiffness of such bearings. The unique 
characteristics of the nanopositioning actuator, i.e., predictable 
and re-configurable open-loop positioning resolution, were also 
demonstrated. Two approaches were presented on how to re- 
configure the open-loop positioning resolution. Advantages and 
disadvantages of each approach were discussed in detail. With the 
FELA, the future work will focus on the implementation of this actu- 
ator on applications such as stem-cell manipulation and fiber optics 
alignment. 
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